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ABSTRACT: Space gravitational wave detectors employing laser interferometry between free- 
flying spacecraft differ in many ways from their laboratory counterparts. Among these differences 
is the fact that, in space, the end-masses will be moving relative to each other. This creates 
a problem by inducing a Doppler shift between the incoming and outgoing frequencies. The 
resulting beat frequency is so high that its phase cannot be read to sufficient accuracy when 
referenced to state-of-the-art space-qualified clocks. This is the problem that is addressed in 
this paper. We introduce a set of time-domain algorithms in which the effects of clock jitter are 
exactly canceled. The method employs the two-color laser approach that has been previously 
proposed, but avoids the singularities that arise in the previous frequency-domain algorithms. In 
addition, several practical aspects of the laser and clock noise cancellation schemes are addressed. 



I. Background 



A. Introduction 

During the last ten years, work has been going on to define and design a spaceborne laser 
interferometer for the purposes of the detection of low frequency (10 _4 Hz — > 1Hz) gravita- 
tional waves [1-3]. These interferometers work by passing laser signals between widely separated 
spacecraft, reading out the relative phases of the signals, and combining signals between dif- 
ferent spacecraft to eliminate laser phase noise and enhance the gravitational wave signal. A 
typical design concept is shown in Figure 1. Three spacecraft (S/C) move on trajectories that 
keep them at the vertices of an equilateral triangle, and signals are passed in both directions 
along each of the three long arms thus formed. The passage of a gravitational wave through the 
system will create small changes in the curvature of space through which the laser signals are 
passing, thereby advancing or retarding the phases of the laser signals. The amplitude of the 
phase shift produced is proportional to the amplitude of the gravitational wave, proportional to 
the distance between the spacecraft, and dependent on the orientation of the arm relative to the 
direction of propagation of the wave. The detection of the phase differences between the signals 
in the arms constitutes the detection of the gravitational wave, and the detailed waveform ob- 
served for the phase change reveals information about the astronomical source that created the 
wave - thus providing gravitational astronomy observations of the accelerated massive bodies 
that emit gravitational waves. 

The difficulty in gravitational wave detection is the smallness of the phase change expected 
from the reasonable astronomical sources one might expect to see. In one mission concept, the 
Laser Interferometer Space Antenna (LISA), the lengths of the arms are 5 x 10 9 m. Despite 
this long baseline, the amplitudes of the expected sources are so small that a sensitivity 1 of 
lOpmHz -1 / 2 , or 10 /xcycleHz _1 / 2 of a 1 fim wavelength laser, has been set as a requirement for 
the detector. Unfortunately, the most stable lasers that can be built have phase fluctuations 
many orders of magnitude larger than this. In ground-based gravitational wave detectors, where 
the phase requirements are more stringent still, the problem of laser phase fluctuations has been 
solved by creating an equal-arm interferometer. In this type of instrument, a single laser signal 
goes down two arms, bounces off the end masses, and is recombined by allowing the returning 
signals to interfere with each other. By maintaining the armlengths strictly equal, the laser phase 
noise cancels when the returned signals are combined. To see how this happens, we consider 
two arms radiating from S/C 1 in Fig. 1. Let the phase of the laser in S/C 1 be <f>i{t) and 
the gravitational wave signal in the arm of the detector between S/C 1 and S/C 2 be h 12 (t). 
Then the signal received at S/C 1, assuming that the signal is sent from 1, transponded with 
no change of phase at S/C 2, and then beat again against the laser in 1, is 

V21 (*) = Mt~ 2Li 2 ) - 0i (t) + h 12 (t) + n 21 (t) (1) 

where n 2 \(t) is the total non-laser phase noise in S/C 1 when tracking S/C 2 and where the 

1 It has long been the practice of the gravitational wave community to characterize noise or 
sensitivity in terms of the "root spectral density". The relation between the variance and the 
spectral density of a time series n is given by (n 2 ) = J A , S n (f)df, where A/ is the bandwidth and 

S n is the spectral density with units (n-units) 2 Hz _1 . The root spectral density is simply y 7 S n (f) 
and has units (n- units) Hz -1 / 2 . Over a bandwidth where the spectrum is flat, the relationship 
between the rms amplitude of n and the root spectral density of n is n rms = y 7 S n (f)y/ Af . 



light travel time in the arm is L\ 2 - The signal received from S/C 3 is written by letting 2^3 
in Eq. 1. When the signals received at S/C 1 from the spacecraft at the ends of the two arms 
are subtracted, the resulting interferometer signal is 

zi(t) = 2/2i (t) - l/3i (*) = Mt~ 2£i 2 ) - Mt - 2L 13 ) + fc 12 (t) - fc 13 (t) + n 21 (t) - n 31 (t) (2) 

If the two arms are exactly the same length (L12 = £13), the laser phase noise will disap- 
pear. However, when the ends of the arms are free-flying spacecraft, the armlengths cannot be 
controlled to maintain L\ 2 = L13. In this case another method must be found. 

B. Unequal- Arm Interferometer Algorithms 

Such a method was discovered by Jim Faller and refined and published in 1996 [4]. In this 
method it was recognized that, since the laser phase noise in 4>\{t) is many orders of magnitude 
larger than any of the other noise sources, the signal in one arm of the interferometer (Eq. 1) 
can be used to determine the laser phase noise, and the measured time series of the noise can 
then be used to correct the interferometer signal (Eq. 2) for the fact that the arm lengths are 
not equal. Working in the frequency domain, one writes the connection between §\ and y 2 \ as 

y2i(f) = (i-e- 2m f L ^)Mf), (3) 

where the Fourier decompositions are given by 

POO POO 

V2i(f)= y2i(t)e 2m f t dt and Mf)= 4>x(t)e 2 ^ 1 dt (4) 
Jo Jo 

The Fourier composition of the laser phase noise is then found by dividing the Fourier decom- 
position of the observed j/21 by the transfer function, 

A ( f\ 2/21 (/) 

Mf) = 1 _ e -2«ifL 12 < 5 > 

and the time series for 4>\(t) is generated by Fourier synthesis from <p\{f). Once 4>\(t) is known, 
its contribution to Eq. 2 can be synthesized and subtracted from the observed z\ (t) to give an 
interferometer signal that is free of laser phase noise. 

The problem in this procedure, however, is that Eq. 5 has singularities at frequencies 
/ = n/L, where n is a positive integer, 2 so the gravitational wave detector will be relatively 
insensitive near these frequencies. In addition, this method requires a Fourier transform, with 
all its sensitivity to biases and aliasing, before the laser phase correction can be implemented. 
Fortunately, a new method has recently been discovered that works entirely in the time domain 
[5] and avoids these difficulties. This method consists of combining the signals from each arm 
with a time-offset in such a way as to undo the inherent time-offsets at the light-times (2Ljj) in 
Eq. 2. The combination is named X{t) and is given by 

X{t) = y 21 (t - 2L 13 ) - y 3 i (* " 2Li 2 ) - y 2 i(t) + 2/31 (t) (6) 

2 Eq. 5 is also singular at / = (when n = 0) , but the laser phase noise in the interferometer 
signal (Eq. 2) goes to zero as / — > 0, so the detector sensitivity remains unimpaired in the low 
frequency limit. 



A little algebra will show that the combination of signals in Eq. 6, using Eq. 1 for Uij(t), 
will exactly eliminate the <f>i(t) terms. The in Eq. 6 represent round-trip signals. In a 
recent paper, Armstrong, Estabrook and Tinto [6], have used one-way signals and identified 
several more combinations of signals that will likewise cancel out the (f>i (t) noise terms without 
cancelling the gravitational wave signals. 

C. High Doppler Rate Algorithms 

As a result of these procedures, laser phase noise can be essentially eliminated as a noise 
source. However, there remains another noise source, whose level is likewise many orders of 
magnitude greater than the desired noise floor, that must be addressed. This noise arises due 
to the fact that the relative velocity of free-flying spacecraft will produce a Doppler shift in 
the frequency of the signal received at each spacecraft, so that the beat frequency between the 
received laser signal and the local laser signal will amount to tens of MHz. In addition, if the 
lasers in the two spacecraft are each independently stabilized by their own Fabry-Perot cavities, 
then the frequencies of the two lasers, being determined by the lengths of the cavities which 
cannot be made exactly equal, will differ by even more, probably by several hundred MHz. The 
gravitational wave will appear as a tiny (~ 10/UcycleHz -1 / 2 ) shift in the phase of this radio- 
frequency (RF) beat signal. As we discuss in the next section, the measurement of phase to 
this precision in a signal at this high a frequency will require a frequency standard with relative 
frequency stability Av/v = 10 -17 over a time equal to the period of the gravitational wave we 
are trying to detect (~ 1000 s). This is beyond the capability of the best laboratory frequency 
standards, or 'clocks', and well beyond the capability of those that can reasonably be used in 
space. So what is to be done? 

In a companion paper [7] to the frequency-domain unequal-arm algorithm paper [4], a 
frequency-domain procedure was presented that would eliminate this clock noise. This method 
required that a second laser signal be used along the arms of the interferometer. Each spacecraft 
must therefore broadcast two laser signals, the main signal and a second signal that is frequency- 
offset relative to the main signal by an amount tied to the local clock. Thus, whatever phase 
reference is being used at one spacecraft is sent off to the far spacecraft where it is measured 
and recorded. Then, working in the frequency domain, an inversion like that in Eq. 5 may be 
performed, and a time series giving the jitter in each clock may be formed by Fourier synthesis. 
In this way, whatever error is induced by reading laser phase relative to a noisy frequency 
standard may be simulated and corrected. This frequency-domain method, however, suffers all 
of the limitations noted above for synthesizing 4>i(t), including the poles at / = n/L. 

It is the main purpose of this paper to present a set of time-domain algorithms that may be 
used to process the dual-frequency laser data and eliminate the noise due to clock jitter. These 
algorithms will be developed in Section II and will represent the clock-jitter counterparts of the 
laser phase noise algorithms reported in references [5] and [6]. In Section III, we also present 
several practical aspects of the data processing that have not previously been addressed. 



II. Time-Domain Algorithms 



A. Instrument Concept 

We assume an instrument labeled as in Fig. 1, with three identical spacecraft numbered 1, 
2, and 3. The signal received at time t by S/C 1 from S/C 2 will be denoted as U2i(t)- In the 
previous section, our expressions for the signals assumed that the lasers in S/C 2 and S/C 3 were 
phase- locked to the incoming signals from S/C 1, so that whatever phase S/C 2 received was 
simply bounced back toward S/C 1. While this assumption simplifies the formulas in the case 
of a single interferometer with a single central spacecraft, there are reasons [6,8] for keeping the 
lasers in each spacecraft independent of the others. Indeed, if data from a second interferometer 
(with vertex at one of the other spacecraft) are to be collected simultaneously with the first, 
then the forms of the signals for the second interferometer are not the simple expressions given 
in Section I. In the remainder of this paper we will assume that the laser in each spacecraft is 
locked only to its own Fabry-Perot cavity. The case where some lasers are locked to incoming 
signals may be recovered as a special case of the expressions we derive (see Section III.D). 

We write the phase of the signal sent by S/C i and received by S/C j as 

Vij{t) =<t>i{t- Lij) ~ <t>j{t) + hij(t) + nij(t) (7) 

The notation is the same as in Eq. 1, except that hij and have changed their meaning slightly. 
In Eq. 1, these quantities were the total signal and the total non-laser noise generated during the 
round-trip of the signals, including whatever noise S/C % contributed while it was transponding 
the signal. Here, in Eq. 7, represents the gravitational wave signal generated during the 
one-way trip from S/C i to j, and represents the noise that arises one-way, dominated at low 
Fourier frequencies by equal amounts of position noise in both spacecraft and at high Fourier 
frequencies by shot noise in the laser receiver at S/C j [9] . We further separate the laser phase 
into a part that is a pure constant frequency and a part that is random phase noise 

<f>i{t) = Vit + Pi (t) (8) 

where Vi is the constant laser frequency in S/C i and Pi(t) is the phase noise in this laser. If we 
also write the armlength as an initial distance plus a constant velocity, 

Lij = Lijfi + Vijt (9) 

then Eq. 7, with Eqs. 8 and 9 included, becomes 

Vij{t) = [{vi - Vj - VijUi] t + pi(t - Lij) - pj(t) + hij(t) + nij(t) (10) 

where we have dropped the constant phase offset associated with 1^,0- 

Now let us estimate the sizes of the various terms in Eq. 10. The magnitude of the first term, 
Vi — Vj, will depend on how closely tuned the frequencies of the lasers in the two spacecraft may 
be, considering that each will be locked to its own Fabry-Perot cavity. For the OMEGA mission 
[3], which proposed using independent Nd:YAG (1064 /im wavelength) lasers in each spacecraft 
in this way, a reasonable estimate was found to be 300 MHz. In any event, the next term, VijVi, 
will be of order ~ 10 MHz, so the receiver design must accomodate this fundamental frequency 



for yij, even in the case of non-independent end lasers. The size of the pi terms in Eq. 10 will 
depend on how accurately the lasers are locked to the cavities in each spacecraft (see Section 
III. A for an estimate). The smallest detectable hij(i) gravitational wave phase signal is required 
to be 10 //cycles Hz -1 / 2 , so uncancelled Pi(t) laser noise and other riij(t) phase noise must be 
less than this. Therefore, in order to detect a gravitational wave with period r = 1000 s in a 
bandwidth Af = 1/r, a phase shift of 5(f) = y/~S$\/Af = 0.3/icycle must be detected in a total 
phase of 300 MHz x 1000 s = 3 x 10 11 cycles. If we recast this requirement in terms of frequency, 
we see that this corresponds 3 to detecting a 2 nHz frequency shift in a signal of 300 MHz. In 
order to measure a frequency to a this accuracy, a frequency standard with frequency stability 
2nHz/300MHz « 10" 17 is required. Since there presently exists no space-qualified frequency 
standard with such performance, clock jitter will be a major noise source in the detectors unless 
it is dealt with in some other way. 



B. Clock Jitter 

To see how clock jitter noise enters the readout of the laser phases, let us consider the 
detection process in more detail. We will actually so far anticipate our final solution to this 
problem as to discuss the two laser signals that will ultimately be needed. Figure 2 shows 
symbolically the signals that are sent and received by the two spacecraft at the two ends of a 
single arm. The fundamental laser frequency in S/C 2 is v 2 and the laser phase noise is p 2 (t) . 
In addition, a second laser signal is superimposed on the same beam, either by merging signals 
from two separate lasers (as in OMEGA [3], with the beat frequency between the two acting as 
the local spacecraft clock) or by modulating the main laser signal at a frequency equal to the 
fundamental frequency of the ultra-stable oscillator (USO) that serves as the local spacecraft 
clock (as in the current LISA design [2]). The frequency of the local clock is f 2 and it is assumed 
to have phase jitter q 2 (t). Thus S/C 2 will send laser signals v 2 t+p 2 {t) and (1*2+ f2)t+P2(t)+q2(t) , 
while S/C 1 has local laser signals V\t + Pi(t) and (ui + fi)t + p\(t) + q\(t). The phases of the 
two signals from S/C 2, as received at S/C 1, are 

4> 2 {t) = u 2 {l-V 12 )t + p 2 {t- L 12 ) + h 21 {t)+n 21 {t) (11a) 

and 

<f)' 2 {t) = {u 2 + / 2 )(1 - V 12 )t+p 2 (t - L 12 ) + q 2 (t - L 12 ) + h 21 (t) + n' 21 (t) (116) 

The n' 21 in Eq. lib may include some noise, like spacecraft position noise, that is the same as 
in n 21 , and other noise, like shot noise, that will be different from the shot noise part of n 21 . 
When these signals are beat against the local laser at S/C 1, there result the two signals 

Ifei(t) = Mt) - Mt) = [iyz ~ v\) ~ V 12 V2] t+p 2 {t - L 12 ) - Pl (t) + h 21 (t) + n 21 (t) (12a) 
and 

V2l(t) = &{t) ~ = [fa - Vl) + (h ~ fl) ~ Vl2{»2 + / 2 )] t 

+ p 2 (t - L12) - pi(t) + q 2 (t - L12) - qi(t) + h 21 (t) + n' 21 (t) 

both composed of constant-frequency terms plus phase fluctuations. 

3 The spectral density of frequency noise is related to the spectral density of phase noise by 
S v = u) 2 Sj, = 47r 2 / 2 S^, so the rms frequency noise is 5v rms = \fS^\/Af = 2irf 5(p ims . 



The procedure for reading out the phases of the two signals y 2 i (t) and y' 21 (t) may best be 
understood by reference to Fig. 3. The two incoming signals from S/C 2, fa and fa 2 , fall on 
a photodiode on board S/C 1, where they interfere with a portion of the two signals, fa and 
fa, that are being broadcast from S/C 1. If, as has been proposed [7], the RF clock signals fi 
are at much higher frequency than the 300 MHz of v 2 — v\, then y 2 \ (t) and y' 21 (t) will be the 
only signals in this frequency band (see reference [7], Eq. 14) and will not be confused with any 
other of the six beat frequencies on the photodiode, generated by mixtures of the four input 
frequencies. 

In order to read out these RF frequencies, a local oscillator (LO) is first used to beat the 
signals down to a low baseband where they may be sampled and fit to determine an average 
phase. For the clock jitter cancellation procedure to work, the LO must be phase-locked to the 
local fi clock, as we show below. The LO could be a piece of hardware (like the phase-stable 
frequency synthesizer developed for the OMEGA mission [10]) or a piece of software in a digital 
processor (as long as the clock cycles for the processor are tied to the spacecraft time standard). 
Whatever the realization of this frequency subtractor, however, it is important for the method 
of this paper that the same frequency be subtracted from both y 21 (t) and J/21 (*)- 4 

The LO will produce a frequency that is close to the frequencies of the two signals 2/21 (*) 
and y 2 i{t). This frequency will be tied to the local spacecraft clock by making the LO frequency 
some rational fraction 021 of the local clock frequency f\. The functional forms of the two 
baseband signals coming out of the mixer will then be 

S21 (t) = y 2 i (t) - 021 [fit + qi (*)] 

= [{v 2 - fi) - V 12 i>2 ~ 021/1] t + p 2 (t - L12) -pi(i) - a 2 iqi(t) + h 21 (t) + n 2 i (i) 

and 

*2l(<) =l/2l(*)-021 [fit + qi(t)} 

= [iy2 ~ vi) + (/ 2 - /1) - V 12 (v 2 + h) ~ a 2 i/i] t (136) 
+ p 2 {t - L12) -Pi(t) + q 2 (t - Li 2 ) - (a 2 i + l)gi(t) + /»2l(*) +"2l(*) 

where 021 has been chosen so that 021/1 ~ ^2 — ^1 — ^12^2- Eq. 13a is the main signal, which will 
be present even if there is no second laser frequency in the beam. The clock jitter problem we 
are trying to solve is apparent in this equation. When any laser signal is detected and read out, 
the qi(t) jitter in the local clock in the receiving spacecraft will produce a phase noise 021*71 (i) 
in the result. Assuming a standard space-qualified USO, this term will correspond to frequency 
noise ~ 1 fiHz. This is three orders of magnitude larger than the 2nHz accuracy needed for the 
gravitational wave requirement. Therefore, even after the X(t) combination is formed and the 
Pi(t) terms drop out, there will be aijqi(t) terms that will remain and will dominate the noise 
in the detector. 



4 An example of a receiver that does not satisfy this requirement would be one that processed 
the RF signals directly to measure their phase. Such an instrument can be thought of as an LO 
that beats each signal separately down to DC. In order to do this for both y 21 (t) and y' 21 (t), 
a different frequency would have to be mixed with j/21 (*) than with y' 21 (t), and the clock jitter 
cancellation procedure we describe here will not work at the required accuracy. 



C. Laser and Clock Noise Cancellation Algorithms 



Before we discuss the algorithms that will cancel the qi(t) clock jitter noise terms in Eq. 
13, let us review the combinations of the y%j(t) signals that are used to cancel laser phase noise. 
We summarize the results from the papers by Armstrong, Estabrook and Tinto [5,6], but we use 
notation (for numbering spacecraft and arms) consistent with Fig. 1 and Eq. 7. 

Consider, then, the combination of Dij{t) defined by 

X{t) = yi 2 (t - L w - 2L 13 ) - yi 3 (t - L 13 - 2L 12 ) + y 21 (t - 2L 13 ) 

- 2/31 (* - 2L i2 ) + 2/13 (* - £13) - Vl2(t - L12) + 2/31 (*) - 2/21 (*) 

It may be noted that X(t) is formed from signals that travel only along the two arms L12 and 
L13, with no data included from the L23 arm. If there were ever a failure in the spacecraft such 
that one arm of the interferometer would be no longer operative, that arm could be designated 
the L23 arm, and the other two arms could still be used to provide a laser-phase-noise-free X(t) 
signal. Functions analogous to X(t) may also be formed by permuting each subscript in Eq. 14 
(1 —>■ 2 — > 3) to define a signal Y(t) which requires no data from L 13 , 

Y{t) = y 23 (t - L 23 - 2L 12 ) - 2/21 (t - L12 - 2L 23 ) + y 32 (t - 2L 12 ) 

(146) 

- 2/12 (* - 2L23) + 2/21 {t - L 12 ) - y 23 (t - L 23 ) + 2/12(0 - 2/32W 
and again to give Z(t) which requires nothing from L12. 

Z(t) = 1/31 (t - L13 - 2L 23 ) - 2/ 32 (t - ^23 - 2L 13 ) + y 13 (t - 2L 23 ) 

- y 23 (t - 2L 13 ) + y 32 (t - L 23 ) - y 31 (t - L 13 ) + y 23 (t) - y 13 (t) 

To see how the these combinations work to eliminate laser phase noise, we expand Eq. 14a, using 
Eq. 7 to express the yij{t) in terms of their elements, and find 

X(t) = Mt~ 2^12 - 2L 13 ) - Mt ~ L 12 ~ 2L 13 ) + h 12 (t - L12 - 2L 13 ) + n 12 (t - L 12 - 2L 13 ) 
-Mt~ 2^13 - 2Li 2 ) + ^ 3 (t - L13 - 2L 12 ) - fei 3 (/ - L 13 - 2L 12 ) - n 13 {t - L 13 - 2L 12 ) 
+ 4> 2 {t - 2L 13 - L12) -Mt~ 2Li 3 ) + /n 2 (t - 2L 13 ) + n 21 {t - 2L 13 ) 

- 03 (t - 2L12 - L13) +<fn(t- 2L l2 ) - h 13 (t - 2L l2 ) - n 31 (t - 2L 12 ) 
+ Mt~ 2Li 3 ) - <fa(t - L 13 ) + h 13 (t - L13) + n 13 (t - L13) 
-Mt~ 2L i2 ) + <fo(t - L12) - h 12 (t - L12) - n 12 (t - L12) 

+ fa(t- L13) - 01 (t) + h 13 (t) + n 3 i(t) 

- Mt - L 12 ) + 01 (t) - ^i 2 (i) - n 2 i(i) 

(15) 

As may be seen by inspecting the first two terms on each line of Eq. 13, the laser phase noise 
terms cancel in pairs. Similar expansions may be used to see the cancellation in Y(t) and Z(t) 
by permuting 1 —> 2 —> 3 in Eq. 15. 

However, as we saw in the last section, it is not the yij(t) signals that are directly measured 
on each spacecraft, but, rather, the Sij(t) signals which contain clock jitter noise in addition to 
the other noise terms. Because one must use the observable Sy(t) to form X(t), there will be 
unavoidable clock jitter noise included from each of the spacecraft. Thus, the creation of 

X(t) = s 12 (t - L 12 - 2L 13 ) - s 13 (t - L13 - 2Li 2 ) + s 21 {t - 2L 13 ) 

(16) 

- s 31 (t - 2L 12 ) + s 13 (t - L 13 ) - s 12 (t - L12) + s 31 (t) - s 21 (t) 



yields a signal with clock jitter 



X(t) = a 12 [qi{t - L 12 ) - q 2 (t - L 12 - 2L 13 )] - a 13 [q 2 {t - L 13 ) - q 3 (t - L 13 - 2L 12 )} 
+«2i [qi(t) -qi(t- 2L 13 )\ - a 31 [q^t) - q x {t - 2L 12 )} + signal + noise 

where the signal and noise terms are the combinations of the hij (t) and riij (t) terms given in Eq. 
15, and where terms representing constant phase shifts have been dropped. The permutations 
1 — ► 2 — > 3 in Eq. 17 give the residual clock jitter noise terms in Y(t) and Z[t). 

The key to the elimination of the a^qk noise terms in Eq. 17 is found by reference to Eqs. 
13. When the s'-it) secondary laser signal (Eq. 13b) is subtracted from the main laser signal 
Sij (t) (Eq. 13a) the resulting difference 5 

r'ij(t) = Sij (t) - 4(t) = (ft - fi + VijStf + q 3 (t) - qi (t - Lij) + n l3 {t) + n'^t) (18) 

will contain a combination of the two clock jitters. The clock frequencies Si wm be- known a 
priori and the orbits will be known well enough to determine Vij. The constant frequency part 
of Eq. 18 may therefore be subtracted off by hand, leaving a signal that contains only the qi 
terms plus the instrumental noise terms: 

rij(t) = r'^t) - (Si - Sj + = qj(t) - Ql (t - L l3 ) + n^t) + < ? (t) (19) 

A little algebra will then verify that the combination 



£(t) = X(t) - a 12 r 21 (t - 2L 13 ) + a 13 r 31 (t - 2L 12 ) - (a 12 + a 2l )r 13 (t - L 13 ) 

+ (a 13 + a 31 )r 12 (t - L l2 ) + (a 12 + a 13 + a 31 )r 21 (t) - (a 13 + a 12 + a 21 )r 31 (t) 

will exactly cancel out the % terms in Eq. 17. Similarly, combinations given by 

ip{t) = Y(t) - a 23 r 32 (t - 2L l2 ) + a 21 r 12 (t - 2L 23 ) - (a 23 + a 32 )r 2l (t - L 12 ) 

+ («i2 + a2i)r 23 (t - L 23 ) + (a 23 + a 12 + a 21 )r 32 (t) - (a 21 + a 23 + a 32 )r l2 (t) 



(20a) 



(206) 



and 



C{t) = Z(t) - a 31 r 13 (t - 2L 23 ) + a 32 r 23 (t - 2L 13 ) - (a 13 + a 31 )r 32 (t - L 23 ) 

+ (023 + a 32 )r 3 i(t - L 13 ) + (a 31 + a 23 + a 32 )r 13 (t) - (a 32 + a 13 + a 31 )r 23 (t) 



(20c) 



will produce signals based on Y(t) and Z(t) that are cleaned of clock jitter noise. It is also 
apparent that, if the a^- coefficients are much less than unity, the additional noise (n^ and n^) 
contributed by adding in the signals will be negligible compared to that already present in 
the signals X(t), Y(t), and Z(t). In fact, the s\- signals only appear multiplied by a^, so, with 
dij small enough, their phase noise n\ ■ can actually be much larger than without adding 
appreciably to the noise in the final signals. 

In order for the phase combinations of Eq. 20 to be calculated, the ratios a^j will have 
to be known. If these ratios are determined via phase-locked loops on board each spacecraft, 



5 It is here that the need for a common LO frequency, a^/j, mixed with both the y^- and 
the ylj signals becomes apparent, for only then will the a^qj terms cancel exactly, leaving the 
simple combination of qj given in Eq. 18. 



then those values will have to be telemetered from each spacecraft so that the combinations in 
Eqs. 20 can be formed. However, since the LO on each spacecraft is only required to beat the 
frequency down to a baseband of, perhaps, a few kHz, accurate enough values for a^- may easily 
be determined on the ground from knowledge of the clock frequencies and of the spacecraft 
orbit. In this case the proper ratios, including Doppler shifts, may simply be uploaded from the 
ground and stored in memory on board each spacecraft. 

In addition to the X(t), Y(t), and Z(t) signals, the recent paper by Armstrong, Estabrook, 
and Tinto [6] has identified several other combinations of the six one-way signals from the three 
arms of the interferometer that are also free of laser phase noise. One such combination 6 is 

A{t) = y 3 i(t)-y 21 (t)+y 23 (t- L 13 )-y 32 (t- L 12 )+y 12 (t- L 13 - L 23 )-y 13 (t- L 12 - L 23 ) (21) 

Other laser-noise- free signals may be found by permuting indices 1 — ► 2 — > 3 to give Bit) and 
then C(t). As before, it is not the y^ that are observable, but the Sij. Using Eq. 13a for the 
Sij, Eq. 21 becomes 

A(t) = a 2 iqi(t) - a 12 q 2 (t - L l3 - L 23 ) + a 13 q 3 (t - L 12 - L 23 ) 

- a 3 \qi{t) + a 32 q 2 (t - L 12 ) - a 23 q 3 (t - L 13 ) + signal + noise 

An important difference between Eq. 22 and Eq. 17 may be noted. In Eq. 17, each multiplies 
a difference between two qj terms taken at different times. This is the property that made 
it possible to find a linear combination of the (Eq. 19) that reproduced the difference and 
corrected for it. In Eq. 22, however, each ajj multiplies a single qj term, and there is no way, 
short of the frequency-domain method of Reference [7], to determine a single qj time series 
by itself. Therefore, there do not exist simple time-domain combinations of signals that will 
eliminate clock jitter from A(t), B(t), and C{t). 

Nevertheless, let us consider the following combination of A(t) and signals: 
a(t) = A(t) + a 12 r 23 (t - L 13 ) - a 13 r 32 (t - L l2 ) + (a 12 + a 23 )r 31 (t) - (a 13 + a 32 )r 21 (t) (23a) 

along with its permuted counterparts 

(3(t) = B(t) + a 23 r 3l {t - L 12 ) - a 21 r 13 (t - L 23 ) + (a 23 + a 31 )r 12 (t) - (a 21 + a 13 )r 32 (t) (236) 

and 

7(i) = C(t) + a 31 r 12 (t - L 23 ) - a 32 r 21 (t - L 13 ) + (a 31 + a l2 )r 23 (t) - {a 32 + a 21 )r 13 (t) (23c) 

When Eqs. 19 and 22 are used in Eq. 23a, we find 

a(t) = — (ai2 - a 2 i + a 23 - a 32 + a 31 - a 13 )qi(t) + signal + noise (24) 

Now let us remember that the are determined by a^- ~ (i/i — Uj — Vij )/fj- When the fj clock 
frequencies are independently set on each spacecraft, there is no reason they should have any 

6 This signal is denoted a(t) in reference [6], but we would like to use Greek letters for the 
signals after they are cleaned of clock-jitter, so we use A(t) for the laser-phase-noise-free signal 
given in Eq. 21 and reserve a(t) to represent its clock-jitter-cancelled counterpart. 



particular relationship to each other, and there is therefore no reason for the combination of 
in the coefficient of qi(t) in Eq. 24 to be particularly small. On the other hand, a careful choice 
of these frequencies might accomplish just that. One such choice is to make each proportional 
to the laser frequency on that spacecraft. In fact, this condition is automatically satisfied when 
the local clock is formed by beating together two lasers that are phase-locked to nearby modes 
of the same cavity. Whether the proportionality is achieved via two lasers or simply by tuning 
each USO to a frequency = Xvi where A is some constant, the coefficient of qi(t) in Eq. 24 
becomes 

ai2 ~ a 2 \ + a 23 - a 32 + a 31 - a 13 =[(1 - ^3)^1(^2 + ^3) (^3 - ^2) 

+ (1 - V 13 )is 2 (vi + iajXi/i - 1/3) (25) 
+ (1 - V 12 )v 3 (vi + ^2)^2 - fi)]/(Ai/i v 2 v 3 ) 

The terms containing will each be of order V^(z/j — Vj)/ fi ~ Vijdij, which, assuming typical 
values ciij ~ 0.1 and ~ 10~ 7 , will make contributions of order 10~ 8 . The terms without Vij 
add up to (u 2 — vi)(v3 — ^2X^3 — ^i)/(A^i^2^3), which is small of order 10~ 13 . Thus, the a(t) 
given by Eq. 23, while not eliminating the g, terms of Eq. 22 identically, can nevertheless reduce 
them to a level where they are negligibly small. Similar derivations and similar conclusions hold 
for the f3(t) and j(t) given in Eqs. (23b) and (23c). 



III. Discussion and Applications 



A. Time Resolution Requirements 

The algorithms for producing laser phase noise cancellation (Eqs. 14) and those for can- 
celling clock jitter (Eqs. 20) require that the observed signals be available at exactly the right 
times so that they may be properly combined. There are two parts to this requirement. First, 
the light times must be known with sufficient accuracy and, second, the signals must be 
available with sufficient time resolution for the signal with the correct time offset to be found. 
The light times will only be known to some finite measurement accuracy and the and 
signals will only be sampled at some finite time resolution. In this section, we will calculate 
what the requirements for the light-time measurement and signal time resolution will be. The 
following two sections will discuss how the time resolution may be acheived without excessive 
data rate requirements for the spacecraft telemetry systems and how the light time may be 
measured to the required accuracy. 

Let us assume that an error 5L mn is made in our knowledge of the mn th arm of the 
interferometer, or, equivalently, that one of the signals is not known at exactly t — L mn , but only 
at t — L mn — 5L mn . Then a portion of the pj (t) laser phase noise (Eq. 8) will remain uncorrected. 
To see how this arises, we write each of the Pj(t) as a Taylor series, expanding Eq. 13a at time 
t — L mn to read 

s ij{t — Lmn) = Pi{t — Lij — L mn ) +f>i{t — Lij — L mn )8L mn —pj(t — L mn ) —pj(t — L mn )5L mn (26) 

where we have dropped all contributions to Sij except for the pj terms. It is the laser phase 
noise that dominates, so if we have sufficient light-time knowledge and time resolution to satisfy 
the laser phase noise cancellation requirement we will automatically have sufficient for the clock 
jitter cancellation. When Eq. 26 is used for all the terms in Eq. 16, X(t) becomes 

X(t) = (6L 12 - 6L 13 ) Pl {t - 2L 12 - 2L 13 ) + 5L 12 [p 2 (t - L 12 - 2L 13 ) - p 2 {t - L 12 )} 
- 5L 13 \p 3 (t - Li 3 - 2L 12 ) - p 3 (t - L 13 )] 

If the measurements of L\ 2 and L\ 3 are independent and of order 5L, then the relationship 
between the spectral density of laser phase noise in X(t) and the spectral density of pi is 

S x = [2 + 4sin 2 (2wL 13 ) + 4sin 2 (2u;L 12 )] (5L) 2 S Au (28) 

where S<\ u = Sp is the spectrum of frequency fluctuations. The time-domain differences of pj in 
Eq. 27 produced the sin 2 (2u;L) terms as transfer functions in Eq. 28. If we average Eq. 28 over 
all frequencies to, we find 

S x = 7M(5L) 2 S Av (29) 

giving a requirement on AL of 



SL = 0.36 -^=J- = 3.6 x 10 -6 s (30) 
v Sav 

where, in the evaluation of Eq. 30, we have assumed the usual phase noise requirement of 
10~ 5 cycles Hz -1 / 2 and assumed that the spectrum of laser frequency stability observed in the 



laboratory, S Au(f) ~ 1Hz 2 Hz -1 at a frequency of 1Hz, can be continued to lower frequency 
when the Fabry-Perot cavity has the thermal stability it is expected to have in space. If this 
is indeed the case, a time-resolution of 1/is would reduce the residual laser phase noise to a 
negligible level. If this stability is not acheivable in space, a finer time resolution would be 
needed for sampling the Sij and signals, and a more precise knowledge of the light time 
between spacecraft would be required. 



B. Signal Requirements 

As was discussed in the previous section, the signals must be sampled on board each space- 
craft with a 1 /is or better resolution. However, Eqs. 14 and 21 show that the combinations 
needed to eliminate laser noise require signals from all three spacecraft (i.e., the formation of 
X(t) requires si 2 , S13, and S21 — 531)- As we will show in Section III.D, even in the case where 
some of the lasers are locked to the incoming signals, signals from multiple spacecraft are needed 
in order to form the complete unequal-arm interferometer combinations. Thus, the signals read 
out on board each spacecraft will have to be telemetered to the ground or to other spacecraft 
in order to create the required combinations. The straightforward communication of the time 
series for s^- and with [is time resolution, and with enough bits to define the phase to the 
required accuracy, would require a data rate on the order of 100 Mb/s. This is prohibitively high 
for the space missions that have been proposed. Fortunately, it is also unnecessary, as we shall 
now show. 

In order to detect gravitational waves in the LF frequency band, a sample time of At = 1 s 
is all that is required in the X(t) and other signals. The only reason for using a higher sample 
rate than this would be to provide the higher resolution necessary to correctly create X{t) and 
the other signals by combining data from the various spacecraft with the proper time offsets. 
However, one-second samples of X(t) can be generated by averaging a higher resolution time 
series. If we average Eq. 16 over the desired sample time, At, we find: 

X n X(t)dt 

I ( rt„+At rt n +At 

= At I / Sl2 ^ ~ Ll2 ~ 2Ll ^ dt ~ / Sl3 ^ " Ll3 " 2L ^ dt 

{Jin J t n 

rt n +At pt„+At pt n +At { > 

+ / s 21 (t - 2L 13 )dt - / s 31 (t - 2L 12 )dt + / s 13 (t-L 13 )dt 

J t n J tn J tn 

rtn+At rtn+At rtn+At "I 

- / S l2 {t - L 12 )dt + I S 31 {t)dt- I S 21 {t)dt\ 

J tn •> tn •> tn J 

where X n is the n th sample, taken at time t n = nAt. Therefore, if the signals are available 
on each spacecraft with /xs resolution, all that is required in order to allow the X(t) combination 
to be formed is the integral of each of the eight terms on the right-hand side over the At interval. 
In fact, all that is required, for example, from S/C 2 is the combination 



-/ 

At A 



tn+At 

[s 12 (t - L12 - 2L 13 ) - s 12 (t - L 12 )} dt (32) 



Consider, for example, a case with light times L i2 = 15.7763391s and L i3 = 15.8227142 s. At a 
time ^7i84 = 7184 s, S/C 2 would need to have available the average of si 2 from t = 7136.578033 s 
(i.e., t n — L12 — 2L13 with 1 ^s resolution) to 7137.578033s [t n — L i2 — 2Li 3 + At). The average 
of S12 from 7168.223661 (t„ - L 12 ) to 7169.223661 (t n - L 12 + At) would also be needed. The 
difference between these two 1-second averages is the signal required from S/C 2 in order to 
form X(t). 

If we consider the total data requirements from all spacecraft and for the complete set of 
variables, x(t), ip(t), and £(£), the following data requirements table may be formed: 
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Table 1. Signal requirements for the set x(t), ip(t), and ((t). The times listed for each signal 
are those at which the signal must be accumulated on board each spacecraft. 



As may be seen, each spacecraft is required to produce four one-second averages of each 
incoming Sjj signal and three of each signal, each with its particular time offset. Since the 
Tij signals are required at the same times as their respective Sij signals, all that is required to 
produce them is to generate the s'^ at the three times given above for the and to simply 
subtract to get = — s'^ as in Eq. 18. The fourteen pieces of data at each spacecraft 



are then combined into three signals that are telemetered to other spacecraft or to the ground 
in order for the %(*)> ^(*)> an d ((t) set of signals to be formed. The data required from each 
spacecraft are 



S/Cl : 



D lx = s 2 i(t - 2L 13 ) - s 2 i(t) - s 3 i(* " 2L 12 ) + s 3 i(*) 

-a2iJ"2i(t - 2L i3 ) + (a 12 + a 13 + a 31 )r 2 i(t) 
= s 2 i(t - L 12 ) - s 2 i(t - L 12 - 2L 23 ) - (a 23 + a 32 )r 21 (t 



L 



Die = s 3i(t ~ L13) - s 31 (t - L13 - 2L 23 ) - (a 23 + a 32 )r 3 i(t - L 



12) 
13) 



(33a) 



S/C2 : 



( D24, = s 32 (t - 2L 12 ) - s 32 (t) - s 12 (t - 2L 23 ) + s 12 (t) 
-a 3 2r 3 2(t - 2L l2 ) + (a 23 + a 12 + a 2 i)r 32 (t) 
D 2 ( = s 3 2(t - L23) - s 32 (t - L 23 - 2Za 3 ) - (ai 3 + a 3 i)r 32 (i - L 23 ) 
D 2x = s 12 (t - L12) - s 12 (t - L12 - 2L 13 ) - (ai 3 + a 31 )r 12 (t - L 12 ) 



(336) 



D 



'3C 



S/C3 : 



D 



3x 



s 13 {t - 2L 23 ) - s 13 (t) - s 23 {t - 2L 13 ) + s 23 (t) 

-ai 3 r 13 (t - 2L 23 ) + (031 + 023 + a 32 )ri 3 (t) 
s 13 (t - L 13 ) - s 13 (t - L 13 - 2L 12 ) - (a 12 + a 21 )r 13 (t 



Lis) 



(33c) 



D 3 ^ = s 23 (t - L 23 ) - s 23 (t - L 23 - 2L 12 ) - (a 12 + a 21 )r 23 (t - L 23 ) 



It should be noted that all six values of the a^- must be known by each spacecraft in order to 
form these combinations on board. It can easily be verified that x(0j ${P)i an d C(0 are formed 
by adding together just these nine pieces of data, so this set of signals represents the data rate 
requirement for this type of mission. 

The set of data needed to form a(t), (3(t), and 7(2) is summarized in Table 2. 
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Table 2. Signal requirements for the set a(t), /3(t), and j(t). The times listed for each signal 
are those at which the signal must be accumulated on board each spacecraft. 

In the case of the a(t), (3(t), and j(t) combinations, there are fewer time offset signals that must 
be accumulated on each spacecraft (three Sij instead of four and two instead of three), but 
the times are not all the same as those in Table 1. Thus, if both sets, x(t), ip{t), and ((t) and 
a(t), /3(t), and are desired, then more data must be accumulated in the laser receivers (six 
and four r^). The data combinations needed to form the a(t), (3(t), j(t) set are 



S/Cl : 



S/C2 : 



D la = s 31 (t) - s 2 i(t) + (a 2 i +a 23 )r 31 (t) - (a 13 + a 32 )r 21 (t) 
Dip = s 3 i(t - L12) - s 2 i(t - L13 - L 23 ) + a 23 r 3 i(t - L 12 ) 
Dij = s 3 i(t - L12 - L 23 ) - s 2 i(t - L13) - a 32 r 2 i(t - L 13 ) 

D 2/3 = si 2 (t) - s 32 (t) + (a 32 + a 31 )r 12 (i) - (ai 3 + a 2 i)r 32 (t) 
D 2l = si 2 (t - L 23 ) - s 32 (t - L13 - L12) + a 3 iri 2 (t - L 23 ) 
D 2a = si 2 (t - L13 - L 23 ) - s 32 (t - L12) - ai 3 r 32 (t - L12) 



(34a) 



(346) 



S/C3 : 



' D 3l = s 23 (t) - si 3 (t) + (a 3 i + ai 2 )r 23 (t) - (a 32 + a 2 i)ri 3 (t) 
< D 3a = s 23 (t - L13) - si 3 (t - Li 2 - L 23 ) + ai 2 r 23 (t - L 13 ) 
, D 3 p = s 23 (t - L12 - L13) - si 3 (t - L 23 ) - a 2 iri 3 (t - L 23 ) 



(34c) 



As was the case for x(£), ip(t), and ((t), each spacecraft must know all of the a^- from the 
other spacecraft before the ot(t), (3(t), ^{t) combinations can be formed. The data requirement 
for a(t), (3(t), 7(t) is again nine pieces of data, three from each of three spacecraft. The signal 
combinations in Eq. 34 are independent of those in Eq. 33, so, if both the x(*)> ^(*)> an d C(*) an d 
the a(t), (3(t), j(t) data sets are desired simultaneously, the data rata requirement is doubled. 

Assuming that each data point requires, say, 64 bits to accurately represent the phase 
in microcycles, the data required from the spacecraft has been reduced, as a result of the 
considerations in this section, from lOOMb/s to 9 x 64 = 576 b/s. The laser phases will be 
sampled with microsecond accuracy and then averaged to form the fourteen one-second phase 
measurements given in Tables 1 and 2. These are then added together in the proper ratios to 
produce the signal combinations of Eq. 33 or Eq. 34, and this is all that must be telemetered 
to the ground. There remains, however, one important use for the high rate data. This is the 
subject of the next section. 



1 f T ' 2 

C 12 (r) = - / s 21 (t-r)s 12 (t)dt (35) 

1 J-T/2 



C. Measuring the L^- 

Let us consider the cross-correlation function of s 21 and s\ 2 over a record of length T. 

rT/2 
-T/2 

When the Sij are expanded, keeping only the dominant laser phase noise terms, Eq. 35 becomes 

rT/2 

~ T / 2 /•>, 
T/2 , r T/2 {6K>) 



i r 1/z 

C 12 (t) = - / M*-^i2-r)-pi(t-r)][pi(t-Li2)-p2(t)]dt 

1 J-T/2 

I [T/2 1 ,T/2 

= -— / p 2 (t - L 12 - r)p 2 (t)dt - — pi(t - r)pi(t - L 12 )dt 

1 J-T/2 1 J-T/2 



-T/2 1 J-T/2 

where the the p\p 2 terms have been dropped since the statistical independence of p\ and p 2 will 
make then negligibly small. The two remaining integrals in Eq. 36 are just the autocorrelation 
functions of p\ and p 2 with time offsets, 

Ci 2 (r) = -A p2 (t + L 12 ) - A pl (r - L 12 ) (37) 

Thus, the cross-correlation function will have peaks at +L\ 2 and — Li 2 , corresponding to the 
zero-lag peaks in the autocorrelation functions of p\ andp 2 - Examination of the cross-correlation 
function and identification of these peaks will therefore provide a measure of the light time L\ 2 . 

The measurement of the will not need to be performed on a continuing basis, since the 
observed Doppler rate in s^j will give the change in range between spacecraft. The measure- 
ment will therefore only be needed once initially and again whenever a check of the integrated 
range is desired. For each range measurement, the Sij may be sampled at a raw rate of one 
every microsecond and the raw data from several seconds of measurement of and Sji can be 
telemetered to the ground at a low data rate over a long period of time. Once both signals are 
available, the identification of Lij will consist simply in forming Cij and finding the unique value 
of r where the two maxima of the cross-correlation function occur. If the data are averaged over 
several sample times, in order to reduce the data required for forming Cij, then interpolation 
between the two maxima on either side of r = ±£ij can determine the correct value. 



One practical problem in the procedure we have just described is that the solid-state lasers 
we have assumed for the space missions are, in fact, very low in phase noise at a sample time 
of 1 fis. The products of the pi laser phase noises risk being swallowed up in the shot 
noises in the phase readout system. There are two things that may be done to get around this 
problem. First, since there is always a phase modulator in the main laser beam that is used 
for the phase locking to the Fabry-Perot cavity, that phase modulator could be used to impose 
a 1 MHz pseudo-random phase noise code on the laser signal. The resulting phase shifts could 
be as large as desired and could be much larger than the shot noise floor. Alternatively, the 
process of averaging the raw data rate over several samples will give a laser phase noise at what 
is effectively a longer sample time. The noise at this sample time will be greater than that at 
the raw rate of 1 fis since lasers typically exhibit a 1// noise spectrum in this frequency range. 
This noise may be great enough by itself to be detectable above the shot noise, and the true 
value of r could be determined by interpolating the correlation function. In Fig. 4, a simulation 
of this general technique is shown for a 1 ms (instead of a 1 fis) time resolution. A total of 
10 5 points of pseudo-random Pi(t) phase noise data were generated with a fundamental time 
resolution of 0.1 ms. The data were then lagged by the proper amounts to represent si 2 and S21 
to 0.1 ms resolution. The Sjj were each averaged over 10-point bins to form a total of 10 seconds 
of 1-ms-resolution signal. A sample of the time series of s 12 is shown in Fig. 4a. Finally, the 
cross-correlation function Ci2(r) was formed (Fig. 4b), and an interpolated estimate of L 12 was 
made, as given in the figure. A similar set of graphs is shown in Fig. 5 for assumed 1/ / noise. In 
each case, with only 10 4 data points in the cross-correlation, the interpolation has determined 
the correct value of L 12 (133.4 ms) to about 0.4% of the 1ms resolution. 



D. Lasers Phase-Locked to the Incoming Signals 

The implementation of the space interferometer that most resembles the laboratory Michel- 
son interferometer is that where two of the spacecraft act like active mirrors, transponding 
without any change of phase whatever signal they receive. For example, S/C 1 in Fig. 1 might 
contain the master laser and the lasers on S/C 2 and S/C 3 would be phase locked to the signal 
they received from S/C 1. As we have shown in Eq. 2 and Eq. 5, the forms of the Michelson 
signal zi(t) and the interferometer signal X(t) are particularly simple in this case and may be 
completely formed with signals that are available on S/C 1. This simple form may also be found, 
beginning with the fundamental formula for one-way arms (Eq. 14a) by setting y\ 3 and y\ 2 to 
zero, reflecting the fact that the phase-locked lasers on S/C 2 and S/C 3 will be locked in such 
a way that this condition is always met. However, when these two quantities are set to zero in 
the formula for Y(t) and Z(t), the results are 



Y(t) = y 23 (t - L 23 - 2L 12 ) - y 21 (t - L 12 - 2L 23 ) + y 32 (t - 2L 12 ) 
+ V2i(t - L12) - y 23 {t - L 23 ) - y 32 (t) 

and 

Z{t) = y 31 (t - L l3 - 2L 23 ) - y 32 (t - L 23 - 2L 13 ) - y 23 (t - 2L 13 ) 
+ y 3 2{t - L 23 ) - y 31 (t - L 13 ) + y 23 (t) 



(38a) 



(386) 



One point to be noticed here is that Y{t) and Z(t) cannot be written by permuting indices in 
Eq. 6, since there is a fundamental difference between an independent central laser sending out 
two beams that are reflected and returned and a central laser that is itself locked to an incoming 
beam whose phase is independent. Thus, even when the end lasers are phase-locked, the only 



way to form combinations beyond X{t) that are laser-phase-noise free is to use the formulas in 
Eq. 38. As is explicitly shown in Eq. 38a, the formation of Y(t) requires a combination of signals 
from S/C 3, y 23 {t- L 23 - 2L 12 ) -y 2 z{t- L 23 ), and from S/C 1, y 2 i(t- L 12 -2 L23 ) -y 21 (t- L 12 ), 
in addition to the signal, y 32 {t — 2L\ 2 ) — y 32 (t), from S/C 2. Different combinations are required, 
again from all three spacecraft, in order to form Z[t). 

When the lasers in S/C 2 and S/C 3 are locked to the incoming phase of the lasers from 
S/C 1, the correct formulas for the generation of x(i), ^(*)> CW> a (t), P(t), and 7(t) are found 
by setting y 13 and y 12 identically to zero at all times in Eqs. 14 and in Eq. 21 (with its permuted 
extensions). To see the data requirements from each spacecraft, the same s\ 3 and s\ 2 variables 
should be set to zero in Eqs. 33 and 34. As may be seen in Eqs. 33, the data rate from each 
spacecraft is the same, with or without the end lasers locked to their incoming signals. Only 
if the clocks on board each end spacecraft were also phase locked to the incoming difference 
between su and s' u (so that r\ 2 = r i3 = 0), would the two signals D 2x and D 3x in Eq. 33 not 
need to be telemetered from the two end spacecraft. And, in the case of a(t), (3(t), and 7(i) 
there can be no reduction in data rate at all. 

To summarize, there is nothing wrong with locking the end lasers in the way that has been 
proposed. Such locking will probably be required if there is a desire to reduce the LO frequency 
to the minimal ~ 10 MHz that is obtained when the only frequency offsets are the unavoidable 
Doppler frequency shifts. However, this approach provides essentially no simplification of the 
data taking nor reduction of the data rate requirements for the mission. 



IV. Conclusions 



As a result of the algorithms we have presented, we find that it is possible to completely 
eliminate the clock jitter noise at all frequencies for the case of the X(t), Y(t), and Z(t) variables 
by use of the data combinations in Eqs. 20. For the case of the A(t), B(t), and C(t) variables, it 
is possible to reduce the clock jitter noise to a negligible level by use of Eqs. 23 and by a judicious 
choice of the clock frequencies on board each spacecraft. We have also investigated the timing 
precision and resolution required for creation of the laser phase noise elimination and for the 
clock jitter elimination (Eq. 30) and have suggested a method for measuring the time-of- flight 
of the laser signals to this accuracy. Finally, we have specified how the data rate requirements 
for the mission can be minimized by use of appropriate simultaneous time-offset averages on 
each spacecraft. The data that must be generated by the laser receivers is given in Table 1 for 
the x(f), ip(t), and ((t) variables and in Table 2 for the a(t), (3(t), and j(t) set. The data that 
are to be communicated to the ground in order to form the laser-noise free and clock-jitter free 
signals are then given by Eqs. 33 for x{t), ^(t), and and by Eqs. 34 for a(t), (3(t), and j{t). 
From these requirements, we derive a minimum data accumulation rate of 576 b/s for the set 
of all three spacecraft. In the case where two lasers are phase- locked to the incoming signals 
from one master laser, the data rate for the x(£), Y>(i), an d CW variables can be reduced to 488 
b/s (if the local spacecraft clocks are likewise phase locked to the master spacecraft clock), but 
there is no reduction possible at all for the a(t), (3(t), and j(t) set of variables. 
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Fig. 1. Geometry and notation for the constellation of 3 spacecraft 




Fig. 2. Geometry of a single one-way link showing the 
laser frequencies generated on each spacecraft. 
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Fig. 3. Details of the laser receiver showing the four signals falling on the 
photodiode (PD) and the two 300 MHz radio frequency signals y 21 an d y\\ 
that are the beat frequencies between the four. The frequency synthesizer 
produces a frequency a 2 Ji that mixes both y 2l and y' 21 down to baseband 
signals s 2l and s 2l , with a 2l ~ (v 2 - V l - V n V 2 )/f v 
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Fig. 4. (a) Segment of the 10 s time series for pseudo-random phase noise 
and (b) the cross-correlation function for s u and s 21 with L 12 = 133.4 ms. 



200 



400 



600 



800 



t (ms) 




x (ms) 



Fig. 5. (a) Segment of the 10 s time series for llf inherent laser phase noise 
and (b) the cross-correlation function for s n and s 21 with L 12 = 133.4 ms. 



